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ABSTRACT 

For the Antennae interacting galaxy pair, we have obtained high quality, 
fully sampled 12 CO J=l-0 and 3-2 maps of the regions surrounding the nuclei 
and the area of overlap between the two galaxies. The maps possess an angular 
resolution of 15" or 1.5 kpc, so far the highest resolution maps available at 
both the J=l-0 and 3-2 transitions. In addition, 12 CO J=2-l data have been 
obtained for the positions of the two nuclei as well as in part of the overlap 
region with 20" angular resolution. The 12 CO J=l-0, 2-1, 3-2 emission all peak 
in an off-nucleus region adjacent to where the two disks overlap. Use of the 
conventional X factor yields ~ 4 x 10 9 M molecular gas mass in the overlap 
region. It is difficult to understand how such a large amount of molecular gas 
can be accumulated in this region given the relatively short lifetime of molecular 
clouds and the limited period of time for this region to form. 

Line emission at 13 CO J=2-l and 3-2 is detected at selected points in 
the two nuclei and the overlap region. Both the 12 CO/ 13 CO J=2-l and 3-2 
integrated intensity ratios are remarkably high in the overlap region. This is 
the first published case in which such high 12 CO/ 13 CO J=2-l and 3-2 ratios 
are found outside a galactic nucleus. Detailed LVG modeling indicates that 
the 12 CO and 13 CO emission originate in different spatial components. The 
12 CO emission may originate within a non-virialized low density gas component 
with a large velocity gradient. Assuming a CO-to-H2 abundance ratio of 10~ 4 , 
the X factor given by the LVG model is an order of magnitude lower than the 
conventional value for molecular clouds in the Milky Way, but it scales inversely 
as the assumed value for this ratio. Accordingly, we suggest the possibility 
that the strong CO emission in the overlap region of the Antennae galaxies is 
associated with increased radiative efficiency, possibly caused by a large velocity 
dispersion within the individual molecular clouds. 
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A comparison of the CO J=3-2 emission with the SCUBA 850 /im continuum 
in the Antennae galaxies shows that the CO line emission contributes globally 
46% of the 850 /im continuum flux and that the ratio of 12 CO J=3-2 to SCUBA 
850 fim flux varies by a factor of two across the system. After correcting for the 
12 CO J=3-2 contamination, the dust emission at 850 /im detected by SCUBA is 
consistent with the thermal emission from a single warm dust component with 
a mass of 1.7 xlO 7 M . 

Subject headings: galaxies: individual (NGC 4038, NGC 4039) — galaxies: 
interactions — galaxies: starburst — dust, extinction — galaxies: ISM — radio 
lines: galaxies 



1. Introduction 

The remarkable antennae galaxy pair NGC 4038/9 is one of the nearest merging 
systems, and as such it has been well studied at essentially all available wavelengths from 
radio to X-ray (e.g. Hummel & van der Hulst 1986; Haas et al. 2000; Hibbard et al. 2001; 
Bushouse, Telesco & Werner 1998; Nikola et al. 1998; Mirabel et al. 1998; Vigroux 
et al. 1996; Fabbiano, Schweizer & Mackie 1997). Its optical morphology, such as the tidal 
tail features, has been successfully reproduced by N-body simulations (e.g. Barnes 1988; 
Mihos, Bothun & Richstone 1993; Dubinski, Mihos & Hernquist 1996). Therefore it is an 
ideal object for a detailed study of the correlation between the interstellar medium (ISM) 
and the star formation activity in mergers. Of particular interest about this archetypical 
early-stage merging system is that the emission of CO, mid-IR (MIR), far-infrared (FIR), 
submillimeter (sub-mm) and radio 6 cm and 20 cm continuum all peak in an off-nucleus 
region which is coincident with the overlap between the two colliding disks. The recent 
release of HST images (Whitmore et al. 1999) have revealed a population of hundreds of 
compact young clusters, whose ages are estimated to be less than 10 Myr, distributed 
along the edge of a dust lane in the overlap region. MIR (Mirabel et al. 1998; Vigroux 
et al. 1996), FIR (Evans, Harper, & Helou 1997; Bushouse et al. 1998) and sub-mm (Haas 
et al. 2000) images of this region indicate even more intense starbursts that are heavily 
obscured by dust and therefore inconspicuous at optical wavelengths. Huge amounts of 
molecular gas have been found in the whole system as well as in the overlap region (Gao 
et al. 2001). Assuming a conventional CO-to-H 2 conversion factor (the X factor), the global 
star formation efficiency Lir/Mh 2 , however, is only ~ 3.7L Q /M , which is similar to that 
of giant molecular clouds (GMCs) in our Galaxy. All these suggest that the bulk of the 
molecular gas in this system is probably still in an initial turning-on stage of the starburst. 
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This provides us with an excellent opportunity to study the impact of galaxy interactions 
and induced starbursts on the molecular gas properties. 

The main purpose of our study of the Antennae galaxies is to determine the physical 
properties of the molecular gas in the starbursting overlap region and the two galactic 
nuclei. This study will permit the density and temperature to be determined for the star 
forming regions and a more reliable estimate of the mass of the molecular gas in these 
regions. 

We describe a detailed excitation analysis using high resolution multi-line CO data. 
This analysis is based on fully sampled maps at a common resolution of 15" (1.5 kpc) 
so that line ratios can be estimated reliably without the frequently employed ad hoc 
assumptions about the source size and brightness distribution. The transitions used are 
CO J=l-0 and J=3~2 with both the 12 CO and 13 CO isotopes represented at J=2-l and 
J=3-2 at selected points, permitting us to examine a wide range of molecular gas excitation 
conditions, the X factor and the molecular gas mass. Thermal dust emission at sub-mm 
continuum wavelengths provides an estimate of the dust mass. Our CO J=3-2 maps permit 
us to correct more precisely the 850 /im map by Haas et al. (2000) for the contribution by 
this transition to arrive at a revised estimate of the dust mass. 



2. 12 CO J=l-0 Observations with the NRO 45m telescope 

The 12 CO J=l-0 observations at 115.27 GHz were conducted in two observing runs on 
1998 March 23-26 and 1999 January 17-18 using the 45m telescope of the Nobeyama Radio 
Observatory (NRO) in Japan. The 2x2 multi-beam SIS receiver S115Q was used, which 
provides relatively fast mapping of the extended CO emission in NGC 4038/9. This receiver 
has four beams separated by 34" and with FWHM beamwidth of 15" at the frequency of 
CO J=l-0. The system temperature during these observations was in the range 350-700 
K. The wide band acousto-optical spectrometer (AOS) was used as the back end which 
provides a total usable bandwidth of 250 MHz and a frequency resolution of 250 KHz. 
At 115GHz, these correspond to a total velocity coverage of 650 kms -1 and a velocity 
resolution of 0.65 kms -1 respectively. 

Our map is centered at the overlap region at R.A.= 12 /l '01 m 54.8 s , DEC= — 18°52'55" 
(J2000) and includes the two nuclear regions. We used a grid size of 8.5" which is a quarter 
of the beam separation of the 2x2 array. Sixteen telescope pointings, each with 10-20 
minutes integration, were required to provide a complete 64-point map with angular extent 
68" x 68", though slightly under sampled. Our final map comprises three such 64-point 
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maps covering most, but not quite the entire CO emitting region of the two galaxies. The 
spectrometer band was centered at Vlsr — 1560 kms -1 . The data were taken by position 
switching in azimuth using a beam throw of 5'. This beam throw is sufficient to position 
the reference beam well off the source at all map positions. The antenna temperatures were 
converted to main beam temperatures (T" m &) using a main beam efficiency r] m j, = 0.46. As a 
standard procedure, telescope pointing was checked every 1.5 hours by observing SiO maser 
lines from nearby evolved stars with the S40 receiver at 43 GHz. The pointing precision 
was typically 3"(rms) at both azimuth and altitude as the wind speed was low during most 
of our observing time. 

The data were reduced using the NRO software NEWSTAR and then output to the 
Astronomical Image Processing Software (AIPS) package for further analysis. Figure 1 
shows the grid spectra of 12 CO J=l-0. A linear baseline was subtracted from each spectrum 
and the spectrum was smoothed to 10.45 kms^ 1 (3.9 MHz). Most of the spectra were 
taken during the first observing run with the best weather conditions (T S ys = 350 - 600 K). 
Those points in the southwest (bottom-right) region were taken during the second run and 
the noise level was higher. Some positions were observed in both runs for comparison. The 
profiles generally agreed with each other and these have been co-added to generate the final 
spectra. 

3. The JCMT Observations 

Most of the 12 CO and 13 CO J=3-2 and 2-1 data were taken using the James Clerk 
Maxwell Telescope (JCMT) during two observing runs during 1998 March 27 - 28 and 
December 20 - 23. Additional data for the 13 CO J=3-2 were obtained in 2000 January. 
Receiver B3 was used to observe the 12 CO and 13 CO J=3-2 transitions at 345.796 GHz and 
330.581 GHz respectively, while receiver A2 was used to observe the 12 CO and 13 CO J=2-l 
transitions at 230.538 GHz and 220.399 GHz respectively. The B3 receiver has dual 
polarized mixers, permitting an increase in sensitivity by a factor of \[2 when observing 
in the dual-channel mode. The system temperatures measured were T sys ~ 500 K for 
receiver B3 and T sys ~ 400 K for receiver A2. The weather conditions were very good 
during the 1998 March run, and values as low as T sys ~ 280 K were obtained for receiver 
B3. The Digital Autocorrelation Spectrometer (DAS) was employed which provides a total 
usable bandwidth of 920 MHz (802 km s" 1 ) for B3 and 760 MHz (980 kms^ 1 ) for A2. The 
frequency resolution was 0.625 MHz, which corresponds to a velocity resolution of 0.545 
kms" 1 at 345 GHz and 0.84 kms -1 at 230 GHz. The beam shape of the telescope is closely 
approximated by a gaussian with HPBW=14" at 345 GHz and 20" at 230GHz. 
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All the data were taken by beam-switching in azimuth using a beam throw of 120". 
Focus and pointing were monitored frequently by observing bright quasars (3C273 or 
3C279) and the pointing (rms) error was found to be ~ 2" for B3 at 345 GHz and ~ 2.7" 
for A2 at 230 GHz at both azimuth and altitude. The spectral line intensity is calibrated in 
terms of TJ, i.e. corrected for atmospheric transmission and losses associated with telescope 
inefficiencies and rearward scatter. Conversion to T mb was done by using the formula 
^Inb = T\/r] mb) where the main beam efficiency r) mb is 0.62 for B3 and 0.69 for A2 according 
to the JCMT users guide. 

The 12 CO J=3-2 map was based on the region mapped at CO J=l-0 with the 
NRO 45m. We first produced a 7 x 7 grid map of the overlap region centered at R.A.= 
12 >>01 m 54.8 s , DEC= -18°52'55" (J2000), then mapped the nuclear regions of NGC 4038 
and NGC 4039 with a 5 x 5 grid and a 3 x 3 grid, respectively. With a spacing of 7", we 
have fully sampled the majority of the emission in this system. The integration time per 
point is typically 2-8 minutes depending on the signal. The spectrometer band was centered 
at V LS r = 1560 kms" 1 . Measurements of the CO J=3-2 spectral line of IRC+10216 at 
Vlsr — km s -1 were used to estimate the calibration uncertainty by comparison with 
standard spectra in the JCMT archives. This uncertainty is ~ 10%. Finally, we repeatedly 
monitored the spectral line emission of the map center throughout our runs and found 
excellent agreement among all the line profile measurements. 

The spectra were reduced using the JCMT software SPECX and a linear baseline was 
removed from each spectrum. Figure 2 presents the observed 12 CO J=3-2 profiles with a 
grid interval of 7". To permit detailed comparison with the 12 CO J=l-0 map, the 12 CO 
J=3-2 data cube was imported into AIPS, re-gridded to 8.5" sampling, and convolved to 
15" spatial resolution. The re-gridded cube of the 12 CO J=3-2 data were then smoothed 
to a velocity resolution of 10.86 km s -1 (12.50 MHz). The resulting spectra are shown 
superposed on the corresponding CO J=l-0 spectra in Figure 3. Since the 12 CO J=3-2 
map is fully sampled, the re-gridded profiles in Figure 3 are similar to the original ones in 
Figure 2. Figure 4(a) shows the integrated intensities with CO J=l-0 contours overlaid on 
the J=3-2 maps in gray scale form, and Figure 4(b) shows a contour and gray scale map of 
the ratio of integrated line intensities. 

Two spectra of 13 CO J=3-2 were obtained — one at the center of the overlap region, 
which is also the center of the 12 CO maps, and the other at the nucleus of NGC 4038. 
The total integration time for the overlap region and NGC 4038 nucleus was 1.5 hrs and 
1.0 hrs, respectively. In order to obtain a reliable estimate of the line intensity ratio of 
12 CO/ 13 CO J=3-2, we also performed a measurement of the 12 CO J=3-2 line emission at 
the same positions. The differences between the 12 CO J=3-2 spectra observed at different 
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times were within the calibration uncertainties observed for B3 at this frequency. Figure 5 
presents the 13 CO J=3-2 spectra together with the re-observed 12 CO J=3-2 spectra from 
the corresponding positions. 

A 5 x 5 grid map with 10" spacing at the overlap region centered at R.A.= 12 h 01 m 54.8 s , 
DEC= — 18°52'55" (J2000) and one point each at the two nuclear regions were observed 
at the 12 CO J=2 1 transition. The integration time per point was 2-5 mins. 13 CO J=2-l 
data were taken at the two nuclei and at three selected points in the overlap region with 
offsets (0, 0), (0, 10") and (0, -10") relative to 12 CO J=2-l map center. The V LSR was 
1560 kms -1 and the integration time was 40 minutes for each point. Figure 6 presents 
the 13 CO J=2-l profiles overlaid on the 12 CO J=2-l profiles at the positions where both 
measurements were made. The 12 CO J=2-l profiles at other positions can be found in Zhu 
(2001). 

4. Comparison between the 12 CO J = 1-0 integrated intensities with 

previously published results 

Our CO map is concentrated mainly on the nuclear regions of NGC 4038 and 
NGC 4039, and on overlap region centered at R.A.= 12 A 01 m 54.8 s , DEC= -18°52'55" 
(J2000). It covers most of the emission from the galaxies, but its spatial coverage is not 
complete. For example, it does not extend to the western spiral arm of NGC 4038, and thus 
does not cover the CO emission associated with the blue arc of star formation in this region. 
The total CO J=l-0 flux in our map is 2370 ±310 Jy kms -1 , about 5 times greater than 
that (430 Jy kms' 1 ) detected by Stanford et al. (1990), and also 2 times greater than the 
total flux reported by Wilson et al. (2000). This value is in excellent agreement with the 
single dish measurements covering the two nuclei and the overlap region conducted by us 
using the NRAO 12m telescope (Zhu 2001), by Aalto et al. (1995) using the 15m Swedish 
European Submillimeter Telescope (SEST) (FWHM ~ 45"), and by Young et al. (1995) 
using the Five College Radio Astronomy Observatory (FCRAO) 14m telescope (FWHM 
~ 50"). Based on their more complete but much lower resolution (FWHM ~ 55") map, 
covering a 200" x 150" region, Gao et al. (2001) report a total CO flux of 3595 ± 195 Jy 
kms -1 , about 50% higher than in our map. We have compared their measurements with 
ours for the region in common and found good agreement between the two data sets. To 
get an upper limit on the flux of CO emission possibly missed in our map, we assume that 
the unobserved emission within a 200" x 150" region has a low brightness temperature 
equivalent to twice the noise level seen in the 45m map, which is a = 5 K km s -1 . The 
conclusion is that the missed flux under these circumstances would be about 2300 Jy km 
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s 1 . This estimate is indeed comparable to the difference between our measured value for 
the integrated flux and that reported by Gao et al. (2001). 

Adopting a standard CO-to-H 2 conversion factor X of 2.8 xlO 20 H 2 cm -2 (K kms -1 ) 
"\ or M(H 2 )=4.78x(L co /K kms' 1 pc 2 ) M (Bloemen et al. 1986; Strong et al. 1988), the 
integrated intensity in our map corresponds to a molecular gas mass (1.1 ±0.14) xlO 10 M . 
The overlap region ( 35" x 50") contains 910 ±59 Jy kms^ 1 or 42% of the total CO flux, 
indicating ~ 4.2 x 10 9 M Q of molecular gas in this 3.4 x 4.9 kpc region, which is 1.9 
times more than that in the OVRO map by Wilson et al. (2000) in the same region. The 
two nuclei of NGC 4038 and NGC 4039 contain ~ 1.8 x 10 9 M and ~ 9.5 x 10 s M , 
respectively. 



5. Kinematics 

Our CO maps at both J=l-0 and J=3-2 transitions cover mainly the overlap region 
and the two nuclear regions. The goals of our kinematic study are to distinguish different 
overlapping gas components and gain insight into the origin of the large amounts of 
molecular gas in the overlap region. A complete study of the velocity field of the entire 
interacting system would require a kinematic model, and is beyond the scope of this paper. 

From the channel maps of the 12 CO J=l-0 and 12 CO J=3-2 data in Figures 7(a) 
and 7(b), we can identify three components in the overlap region corresponding to super 
giant molecular complexes (SGMCs) in the OVRO maps by Wilson et al. (2000). The first 
component CI (indicated by a star), centered at (±5", ±5"), occupies the velocity range V 
= 1350 - 1630 kms± The second component C2 (indicated by an empty triangle) is at 
(0", -10") and V = 1410 - 1650 kms" 1 . The third component C3 (indicated by a cross) is 
at (—5", ±25") and V = 1410 - 1630 kms -1 . CI and C2 appear to be associated with the 
disk of NGC 4039, whereas C3 appears to be associated with the region of overlap between 
the two disks. CI and C2 are equally strong, while C3 is weaker. All three components show 
very little variation in position over a velocity range of 240 kms" 1 . Closer examination of 
the line profiles (Figures 1 and 2) reveals a double-peaked feature in CI, so this component 
can be further divided into two sub-components. These are spatially resolved into two 
SGMCs in the OVRO interferometer map which are named as SGMC1 and SGMC2 by 
Wilson et al. (2000). The component C2 is also resolved into three SGMCs, named SGMC3, 
SGMC4 and SGMC5 in the OVRO map. The molecular gas in the component C3 is more 
diffuse than CI and C2, and there are no named SGMCs associated with it, though it is 
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also detected in the OVRO map. 

The line profiles associated with C3 are striking, especially in the region near the line 
AS = 20", or DEC(J2000) = -18°52'35". This line is coincident with the southern edge of 
the NGC 4038 disk where it begins to overlap with NGC 4039, and hereafter we denote this 
line as "OV-line" . From the overlaid 12 CO J=3-2 and 1-0 line profiles in Figure 3, it can 
be seen that the profiles for the two transitions generally agree quite well with each other 
in the central and southern part of the overlap region with declination offset A5 < 10" 
from the map center. However, in the region near the OV-line, the J=3-2 transition peaks 
are offset from the J=l-0 peaks, implying at least two velocity components with different 
excitation. More remarkably, the stronger CO J=3-2 peak switches from the high velocity 
to the low velocity feature when crossing the OV-line from south to north, while the J=l-0 
profiles do not show such a dramatic change. This indicates a sharp excitation gradient for 
both the high and low velocity components near the OV-line. As mentioned above, this line 
appears to be the border between the two optical disks. Therefore the degree of excitation 
in the two disks may be very different from one another. 

The velocities in the overlap region are all lower (blueshifted) than the velocity of the 
two nuclei, both of which have a similar velocity V = 1640 kms -1 . The almost identical 
velocity for the two nuclei suggests that the relative motion between these two galaxies is 
in the plane perpendicular to the light of sight, or that the two galaxies are at the point of 
turning around from receding to approaching each other. This pattern is consistent with 
that indicated by the measurements of HI (van der Hulst 1979; Hibbard et al. 2001) and 
Ha (Amram et al. 1992). From the kinematics of the tails and HI velocity field of the spiral 
arms, Hibbard et al. (2001) concluded that NGC 4039 is in front of NGC 4038 (closer to the 
observer) and both galaxies are rotating anti-clockwise. Higher spatial resolution data from 
CO J=l-0 (Stanford et al. 1990) and Ha line emission (Amram et al. 1992) show a velocity 
gradient around the nucleus of NGC 4039, which suggests that the differential rotation 
is redshifted southwest and blueshifted northeast of the nucleus. This rotation pattern 
can partly account for the blueshifted velocity we see in the overlap region, but cannot 
explain the complicated spectral profile in Figure 3, especially for the component C3. The 
blueshifted velocity component in C3 extends continuously from the NGC 4038 nucleus 
down to the disk of NGC 4039. This could be material being accreted onto NGC 4039 from 
the bridge connecting the two disks. Similar gas components with the same velocity range 
in the same region are also seen in the VLA HI image (Hibbard et al. 2001). In addition, 
the overlap region is also the place where the northern tail of NGC 4039 connects to the 
inner disk. The infalling or accreting material can cause strong shocks and greatly distort 
the velocity field in this region. Further n-body simulations of the encounter including gas 
hydrodynamics with high spatial resolution are necessary to explain the kinematics in the 
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overlap region. 

6. Ratios of the integrated line intensities 

We define the ratios of integrated main-beam temperature as follows: 

r-2i = J T mh ( 12 COJ = 2 - l)dV/ J T mh ( 12 COJ = 1 - 0)dV 
r-3i = J T mh ( 12 COJ = 3 - 2)dV/ J T mh ( 12 COJ = 1 - 0)dV 
r-32 = J T mh ( 12 COJ = 3 - 2)dV/ J T mh ( 12 COJ = 2 - l)dV 
= J T mh ( 13 COJ = 3 - 2)dV/ J T mh ( 13 COJ =2- l)dV 
.„ = / T mh { 12 COJ = 1 - 0)dV/ J T mh ( 13 COJ = 1 - 0)dV 
R21 = J T mh ( 12 COJ = 2- l)dV/ J T mb ( 13 COJ =2- l)dV 
R32 = J T mh ( 12 COJ = 3 - 2)dV/ J T mh ( 13 COJ = 3 - 2)dV 

6.1. A gradient in the CO J=3-2/J=l-0 ratio 

After convolving the CO J=3-2 data from 14" to 15" to match the resolution of the 
CO J=l-0 data we can obtain reliable line intensity ratios of 12 CO J=3-2/J=l-0 (r 31 ). 
We used only those points with good S/N ratio for both the 12 CO J=l-0 and 12 CO J=3-2 
spectra. Table 1 lists the positions, integrated intensities at CO J=l-0, 3-2 and the r 3i 
ratios, together with the uncertainties in the latter quantity. Figures 4(a) and 4(b) show 
that the CO J=3-2 peak is offset from the CO J=l-0 peak by ~ 3" in different directions 
in the two nuclei (to the north for NGC 4038 and south for NGC 4039), which may possibly 
be due to a pointing error. This would result in an overestimate of the r 3 i ratio north of 
NGC 4038 and south of NGC 4039 nucleus by approximately 10%. Thus the high r 31 ratio 
in these regions should be interpreted with caution. In the overlap region the two map 
peaks appear coincident, indicating that pointing error should be less than 3". We also 
convolved the 12 CO J=3-2 map to the angular resolution of 12 CO J=2-l (20") to obtain 
the r 32 ratio, which should be less sensitive to pointing errors. The distribution of the r 32 
ratio is similar to that of the r 31 , so the gradient in r 31 seen in the overlap region is most 
likely real. 
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The r 31 ratios in the nuclei of NGC 4038 and NGC 4039 are respectively 0.96 ±0.10 and 
0.77 ±0.10, which is higher than the average value 0.64 ±0.10 in the nuclei of nearby spirals 
and starburst galaxies (Devereux et al. 1994; Mauersberger et al. 2000). The gradient in 
r 3 i across the overlap region is striking. It is as low as 0.3 - 0.4 in the northern part of 
the overlap, but changes to 0.8 — 0.9 toward the southern edge of the overlap region. The 
r 31 ratio also drops by a factor of two toward the east and west edge of the overlap region. 
Figures 8(a) and 8(b) show respectively the r 31 contours overlaid on the K-band NIR image 
(Hibbard et al. 2001) and MIR 15 /im image from ISOCAM (Mirabel et al. 1998). The 
NIR image reveals two young clusters in the southern overlap region which are ~ 4 times 
brighter than the clusters in the northern overlap. The bright knot in the southern overlap 
in the ISO MIR image indicates large amounts of hot dust which is likely to be heated by 
the young clusters there. The derived r 31 ratios are spatially correlated with the strong NIR 
and MIR emission, suggesting that high excitation regions are closely associated with the 
warm gas/dust. 

As mentioned in §5, there is also an excitation gradient associated with different 
velocity components near the OV-line DEC(J2000) = — 18 52'35" (see Figure 3). North of 
this line, the low velocity component has a higher r 31 ratio than the high velocity ones. 
This pattern is also seen in the disk of NGC 4038, east of the nucleus. On the other hand, 
south of the OV-line, including the disk and nucleus of NGC 4039, the r 3i ratio is similar 
for both the low and high velocity components. The most striking profiles are in the vicinity 
of the OV-line, where the low velocity component is undetectable in CO J=3-2 (r 31 < 0.1), 
whereas r 3 i is as high as 0.8 — 0.9 in the high velocity component. Since CO J=3-2 is 
generally associated with warmer and denser regions, the high velocity component (with 
high r 3 i ratio) in this region is most likely to originate in molecular clouds in the NGC 4039 
disk, while the low velocity component could originate in the halo where the temperature 
and density may be too low to generate significant CO J=3-2 emission. Considering that 
the low velocity component corresponds to a high blueshifted velocity, and the fact that 
near the OV-line the two optical disks appear to be overlapping, this component could also 
be the gas in the bridge being accreted by NGC 4039. Detailed modeling including both gas 
dynamics and radiative transfer is needed to test this hypothesis and explain the excitation 
gradient over the velocity range. 



6.2. 12 CO/ 13 CO line intensity ratios 



At some selected points in the overlap region and the two nuclei, isotope emission 
CO J=3-2 and J=2-l was detected making it possible to estimate the i? 32 and R 2 i- 
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These ratios should be relatively free from errors in pointing and absolute flux calibration 
since the two lines were observed immediately one after the other with the same receiver. 
Thus the primary contributors to the uncertainty are thermal noise and baseline errors. 
The derived R32 ratio corresponds to 14" angular resolution. If this ratio holds also for the 
20" resolution associated with R21, and assuming that the 13 CO J=3-2 emission follows 
the distribution of 12 CO J=3-2 emission, we can estimate the 13 CO J=3-2 line intensity 
for a 20" beamsize and derive the 13 CO J=3-2/ 3=2-1 ratio, which is 0.77 ± 0.23 for the 
center of the overlap region and 0.60 ± 0.20 for the nucleus of NGC 4038. Table 2 lists the 
positions, line ratios and their associated uncertainties for the five points with 13 CO line 
measurements. The higher resolution 12 CO J=3-2 emission has been convolved to 20" to 
match the measurements at lower resolution. The convolved r 31 ratio is slightly lower than 
that in Table 1 which corresponds to 15" resolution. 

The largest value of the ratio R 2 i = 29 is at offset (0, +10") which is in the overlap 
region, 34" (3.3 kpc) away from the nucleus of NGC 4039. This value is a factor of 2 
higher than the mean value of R21 for normal spiral galaxies and starburst galaxies (Aalto 
et al. 1995). The value of R32 at the overlap region (0,0) central position is also higher by a 
factor of 3 than those found by Wall et al. (1993) in the nuclear regions of nearby starburst 
galaxies. To our knowledge, this is the first published case in which high R 2 i and R32 ratios 
have ever been found outside the nucleus of a galaxy. 

It is well known that interacting galaxies, especially mergers, tend to have higher 
Rio integrated line intensity ratios than normal spiral galaxies (Aalto et al. 1991, 1995, 
1997; Casoli et al. 1991, 1992a, 1992b; Hurt k Turner 1991; Turner & Hurt 1992; Garay, 
Mardones, & Mirabel 1993; Henkel & Mauersberger 1993; Henkel et al. 1998). The origin of 
this so-called "missing 13 CO" problem (Casoli et al. 1992b) is not yet well understood. Aalto 
et al. (1991, 1995, 1997) explained this effect by a relatively low optical depth at the 12 CO 
J=l transition, but many investigators attribute it to the high [ 12 CO/ 13 CO] abundance 
ratio {e.g. Henkel et al. 1998; Taniguchi & Ohyama 1998, Taniguchi, Ohyama & Sanders 
1999). Such abundances can be produced by (1) the inflow of unprocessed disk gas with high 
12 CO / 13 CO abundance ratios, (2) the 12 C enhancement caused by selective nucleosynthesis 
in massive stars (e.g. Henkel & Mauersberger 1993), (3) selective photodissociation of the 
13 CO molecules in the UV-intense starburst environment (e.g. Fuente et al. 1993). 

From starburst chemical evolution models Henkel & Mauersberger (1993) predict that 
large amounts of low 13 C abundance material should exist in the outer disks of spirals. 
However, the extreme weakness of the 13 CO J=l-0 line makes it difficult to detect in the 
outer region of normal spirals. This problem is alleviated by galaxy interaction which can 
enhance the CO emission on the outer disk through cloud collisions and induced starbursts. 
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Our finding of high R 2 i and i? 32 ratios in the overlap region of the Antennae is consistent 
with the gas infall and chemical evolution scenarios. The region with the highest values 
of i?2i is associated with the youngest starburst clusters in the HST image, whose ages 
are estimated to be less than 5 x 10 6 years (Whitmore et al. 1999). These scenarios can 
also naturally explain the gradient of the i? 2 i ratios, i.e. highest in the north and lowest in 
the southern part of the overlap as well as in the two nuclei. The kinematic data (see §5) 
suggest that most of the gas infall/accretion activity is taking place in the northern part of 
the overlap region near the OV-line at DEC(J2000) = -18°52'35" 

However, we note that many mergers with a high R w ratio do not necessarily have a 
high R 2 \ ratio compared to other starburst galaxies (Aalto et al. . 1995). This leads Aalto 
et al. (1997) to argue that the high R w ratio is associated with small, warm but compact 
molecular clouds with low-to-moderate optical depth (r ~ 1) in the 12 CO J=l-0 line. With 
the measurement of both R 2 \ and _R 32 , we can perform an excitation analysis to better 
constrain the 12 CO/ 13 CO abundance ratio. 

7. Excitation analysis 

The variation of the r 31 , R 2 ±, and i? 32 ratios in the Antennae system indicates dramatic 
changes in excitation conditions. These ratios cannot readily be explained under the 
assumption of LTE. To estimate the physical parameters of the molecular gas, we have 
employed a large velocity gradient (LVG) model (e.g. Goldreich & Kwan 1974). In such 
models it is assumed that the systematic motions rather than the local thermal velocities 
dominate the observed linewidths of the molecular clouds. We use this model to fit the 
observed line ratios in Table 2, for different combinations of (7^ , n(H 2 ), A), where Tk is 
the kinetic temperature and A = Xco/(dV/dR), with Xco= [ 12 CO /H 2 ] being the fractional 
abundance of 12 CO with respect to H 2 and dV/dR the velocity gradient. The optimum set 
of parameters is determined by minimizing x 2 . 

An extensive grid of LVG models was searched with a parameter range of Tkm = 5 — 200 
K (AT kin = 2 K), n(H 2 ) = (0.1 - 10 4 ) x 10 3 cm' 3 (Alogn(H 2 ) = 0.5), A = (0.1 - 10 2 ) x 10~ 6 
pcjkms- 1 ]- 1 (AlogA = 0.5), and 77 = [ 12 CO/ 13 CO] = 20 - 100 (Arj = 5). These ranges 
cover all possible conditions found in the GMCs of our Galaxy as well as in external 
galaxies. The results are summarized in Table 3. Columns 2-4 contain the LVG parameters 
noted above, column 5 contains the CO column density, column 6 the optical depth at the 
J=l transition, column 7 contains the brightness dilution factor (defined as the ratio of 
the observed main beam brightness temperature to the radiation temperature given by the 
LVG model), column 8 contains the assumed isotope abundance ratio, and column 9 shows 
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the value of \ associated with the fit. 

7.1. Single component modeling 

The center position of the overlap region (offset (0,0)) and the nucleus of NGC 4038 
have both 13 CO J=3-2 and J=2-l measurements, which enable us to better constrain the 
physical parameters. At the offsets (0, +10") and (0, —10") in the overlap region and the 
nucleus of NGC 4039 where only R 21 was measured, we assume that the i? 32 ratios at these 
positions are similar to that of the (0,0) position, but with an uncertainty of 50%. Table 3 
lists the estimated physical parameters for these regions using a single component model. 

Table 3 shows that a single component model can fit the line ratios in the NGC 4038 
nucleus reasonably well, with Xmin— 1-2- The best fit corresponds to T k = 43 K, n(H 2 ) = 
3.5 x 10 3 cm" 3 , A = 3.2 x 10~ 6 pc [kms" 1 ]" 1 which gives a r 2 i = 1.1 , r 31 = 0.82, R 21 = 
13.7, i?32 = 20 and 13 r32 = 0.51. However, for the overlap region where the R 2 \ and i? 32 
ratios are very high, no acceptable fit can be found with a single component model. The 
closest fit (xmin = 2.2) corresponds to moderate optical depth gas with r ~ 0.9, 1.5, 3.4 at 
the 12 CO J=l-0, J— 2-1, and J=3-2 transitions respectively (Table 3), with an abundance 
ratio rj = 70. This fit yields 13 r 32 = 0.48, considerably lower than the observed value of 
0.77. Thus unless we have overestimated 13 r 32 by 40%, a single component model cannot 
reproduce the observed line intensity ratios. In general the ranges rj = 40 - 200, T k = 28 - 
160 K, n(H 2 ) = (1 - 6) x 10 3 cm" 3 and X co /(dV/dR) = (1 - 10) xlO" 6 pc [kms^ 1 ]- 1 can 
fit the ratios equally well with x < 3. When r\ is less than 40, no fit is possible. This 
indicates that the probability for a larger rj is high in the overlap region. 

Further to this point, it is probable that two components are necessary to model 
the overlap region. A lower limit on gas density can be set by the 13 r 32 ratio. If the 
13 CO emission is from the same component as 12 CO, the 13 CO J=3-2 must be optically thin 
(since i? 32 >> 1) and normally excited by collisions. The critical density of 13 CO J=3-2 
is relatively high (~4x 10 4 cm -3 ) and moderate gas density is needed to excite enough 
13 CO molecules to the J=3 level. The result is that n(H 2 ) must be greater than 1 x 10 4 
cm -3 to satisfy 13 r 32 = 0.8. However, at this density, the optically thick line 12 CO J=3-2 
is nearly thermalized due to radiative trapping and the r 31 ratio should be close to unity. 
The result is that with a r 3i ratio of 0.78 ± 0.08, the 12 CO emission is most likely from a 
gas density lower than 10 4 cm -3 . This situation is also true for other kinetic temperatures 
we searched from 10 K - 200 K. 

Such discrepancies clearly demonstrate that a single component model fails to describe 
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the complicated physical conditions in the ISM of the Antennae. This is not surprising since 
the optical depths of 12 CO and 13 CO are significantly different and they may be probing 
different regimes of the molecular cloud with different physical conditions. In addition, 
the complex line profiles in the overlap region indicate the presence of more than one 
component. 

7.2. A two-component model for the overlap region 

Our analysis in the previous section suggests the presence of at least two components: 
a low density component [n(H 2 ) < 10 4 cm" 3 ] dominating the 12 CO emission and a high 
density component [n(H 2 ) > 10 4 cm -3 ] dominating the 13 CO emission. This is consistent 
with the findings of other investigators from the studies of other galaxies (e.g. Aalto 
et al. 1994, 1995; Wall et al. 1993). As suggested by Aalto et al. (1995), in regions with 
strong tidal force, such as the overlap region of the Antennae and the galactic nucleus 
with strong differential rotation, gas of low density will not survive as clouds, but will 
be torn apart and sheared into filaments. Dense gas, on the other hand, will remain 
in cloudy structures. Alternatively, the two-phase model can also be a combination of 
photon-dominated regions (PDRs) (e.g. Tielens & Hollenbach 1985) on the surface of 
13 CO-emitting dense cores. The low density phase is more diffuse and has a much larger 
filling factor, which makes it the dominant contributor to the 12 CO emission. The dense 
cores have a much smaller filling factor, but may have significant optical depth (i.e., r ~ 1) 
for 13 CO, thus contributing most to the observed 13 CO emission. Evidence for such dense 
molecular gas has been uncovered in a number of studies where 13 CO is found to be 
optically thick (e.g. White 1997). 

A two-component model will readily fit the observed line ratios. These two components 
are not resolved either in velocity or in space, and we assume the same values of rj and A 
for both components. Taking the low density component A as the one that best fits the 
data in the single component model (Table 3), we find that a component B with n(H 2 ) = 
1 x 10 5 cm -3 and T k = 90 K can fit the observed line ratios within the measurement errors. 
This component has a 13 r32 ratio of 1.7. The combination of these two components with 
dilution factors (defined as the ratio of the observed to model brightness temperatures) ] a 
= 5.1 x 1(T 2 and f B = 3.3 x 10~ 4 gives r 21 = 1.06, r 31 =0.80, R 21 = 23, R 32 = 28, 13 r 32 = 
0.6. In Table 4 we list four combinations that can fit the observed data reasonably well. 

The low density component generates most of the 12 CO emission, with ratios R 2 \ > 26 
and -R32 > 25. The r 32 ratio is 0.7-0.8 in this component, but 13 r 32 is smaller than 0.5. 
Thus the high density component requires 13 r 32 > 1, which implies warm and optically thin 
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3 CO emission. It should have a filling factor much smaller than that of the low density 
component, so that it has little effect on the 12 CO emission. The small dilution factor of 
the dense component is consistent with the weak HCN J=l-0 emission in the overlap region 
and the nuclei of NGC 4038 (Gao et al. 2001). Gao et al. (2001) find that the L HCN /L C o is 
lower for the Antennae compared with other starburst galaxies, implying that the fraction 
of high density gas in this system is smaller than other IR luminous mergers. 

It must be noted that, by introducing a two-component model, there are twice as many 
free parameters and thus too few observed line intensity ratios to constrain them. The data 
can be satisfied with very different parameters T = 30 - 130 K, n(H 2 ) = (1 — 8) x 10 3 cm -3 
for the low density component and T=40 - 200 K, n(H 2 ) > 3 x 10 4 cm" 3 for the high 
density one, with the assumption that both components have the same A = 10~ 6 — 10~ 5 
pcjkms" 1 ]" 1 and isotope ratios = 40 - 100. The LVG fit of the line intensity ratios is 
useful only as a rough guide on the average conditions expected in the region covered by 
the telescope beam (~ 2 kpc here), but it shows that the molecular gas that produces most 
of the 12 CO emission must have a low optical depth, whereas the 13 CO emission probably 
originates in a different component with a much smaller filling factor. In addition, the 
abundance ratio 77 is high compared to Galactic molecular clouds. 

Although the single-component model can fit the nuclear region of NGC 4038, the 
results do not rule out the presence of multiple components here as well. 

7.3. Inferences from the r 31 ratio alone 

For most of the region mapped there are measurements of r 2 \ and r 3 i only (i.e. no 
isotope measurements). However, the r 2 \ ratio is close to unity indicating that the 12 CO 
J=2-l emission is thermalized. Thus it cannot provide much of a constraint on the physical 
parameters. However the r 31 ratios are under unity for all regions, and are sensitive to 
excitation conditions. In the nuclear and southern overlap regions, r 31 ~ 0.8 and this value 
sets a lower limit to the kinetic temperature Tf~ > 20/T since no LVG solution is possible 
otherwise. In the northern part of the overlap region where the r 3 i ratio is found to be 
~ 0.4, the lower limit for the gas temperature is ~ 10 K. When the kinetic temperature is 
greater than the lower limit, the r 31 ratio also depends on the gas density. For T k =20-80K, 
a value of r 3i =0.8 indicates n(H 2 ) ~ 10 3 cm~ 3 , while r 3 i=0.4 implies n(H 2 ) ^ 10 3 cm -3 . 
Near the OV-line the low velocity component has r 3 i < 0.1, this extremely low r 3 i ratio 
suggests that either the molecular temperature is very low ( < 10K ) or the gas density is 
very low (<5x 10 2 cm -3 ). 
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7.4. Molecular gas mass 

In §4 we estimated the amount of molecular gas based on the assumption of the 
conventional Galactic CO-to-H 2 conversion factor. Now we can use the parameters derived 
from the LVG model to get a more knowledgeable estimate of the molecular gas mass. 

From Table 3, the total CO column density can be constrained to a narrow range 
e.g. N(CO)/AV = n(H 2 )( -§f^)f ~ (1.5 - 1.9) x 10 15 cm 2 / kms" 1 for the combination of 
two components in the central 20" of the overlap region. Integrating over the 20" area, we 
have M H2 = (1.4 - 1.7) x 10 8 ( 200 t^ s -i )(^j) M o, wh ere X co is the abundance ratio. In 
the same region, use of the Galactic conversion factor yielded 1.8 x 1O 9 M . Similarly for 
the nuclear region (20") of NGC 4038, we have M H2 = 8.5 x 10 7 ( 100 £L-i )(^) M ©- These 
values depend upon the unknown value for X C o, which is generally assumed to be in the 
range 10~ 5 -10~ 4 for starburst galaxies (e.g. Mao et al. 2000, Booth & Aalto 1998). At the 
lower end of this range, the molecular masses are comparable to those derived from the 
conventional X factor, but X C o = 10~ 4 yields molecular masses which are a factor of 11-13 
times lower. 

Unfortunately, few if any reliable measurements of the abundance ratio exist for 
starburst galaxies. Reliance must be placed on observations and theoretical models of 
Galactic GMC's and star forming regions (e.g. OMC-1 and Sgr B2). Typical values of 
Xco f° r such regions are near 10 -4 , and range from 5 x 10~ 5 to 2.7 x 10 -4 (e.g. Blake 
et al. 1987, Farquhar et al. 1997, Hartquist et al. 1998, van Dishoeck 1998, van Dishoeck & 
Blake 1998). Thus it appears unlikely that X C o < 5 x 10 -5 in starburst regions unless the 
metallicity is abnormally low. If X C o ~ 10" 4 , then the values of A ~ 10~ 6 — 10~ 5 in Table 3 
imply the existence of a low optical depth, such that dV/dR ~ 10 — 100 km s~ l pc' 1 for the 
constituent clouds. Though such large velocity gradients cannot be ruled out by observation 
of individual clouds in starburst systems because of inadequate spatial resolution, they seem 
perhaps physically implausible even for disrupting clouds. It is possible that a combination 
of lower Xco an d large velocity dispersion are responsible for the low optical depths. It 
is likely that such clouds are virially unstable, and that the large velocity dispersions in 
the clouds give rise to a higher CO radiative efficiency (and hence lower X factor). If 
Xco ~ 10~ 4 then the corresponding CO-to-H 2 conversion factors X would be (5.1 - 6.4) 
xl0 19 (|g) cm' 2 [Kkms' 1 }' 1 for the central overlap and 2.3 x 10 19 (^) cm~ 2 [Kkms' 1 ]- 1 
for the nucleus of NGC 4038. These values are 5-13 times smaller than the conventional X 
factor 2.8 x 10 20 cm" 2 [K kms~ 1 } _1 (Strong et al. 1988). These are likely to be lower limits, 
and it seems plausible to adopt a value for X by, say, a factor of 5 lower than the Galactic 
value overall, so that the inferred total mass of molecular gas in the Antennae system would 
be reduced to 2 x 10 9 M(T,. 
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8. Comparison with sub-mm continuum 

Sub-mm continuum data were obtained by Haas et al. (2000) using SCUBA on the 
JCMT. The angular resolution is 15" at 850 /xm and 7" at 450 /xm . The former is essentially 
identical to that for our NRO and JCMT line data. Our CO J=3-2 line data thus permit a 
reliable estimate to be made of the contribution to the SCUBA 850 /xm brightness by this 
line which is contained within the SCUBA bandpass filter. Although SCUBA measures 
primarily the continuum emission from dust, this molecular line contribution must be 
removed before we can interpret the continuum emission. This SCUBA equivalent flux for 
a line in the SCUBA band at the radial velocity of the Antennae galaxies was estimated 
using the SCUBA filter bandpass profile and a beamsize of 15" . The conversion factor used 
for 15" SCUBA beam is 0.735 mJy beam -1 (K km/s)" 1 . Figure 9(a) shows the uncorrected 
SCUBA 850 /xm contours overlaid on the CO J=3-2 greyscale map. It may be seen that 
the SCUBA image closely follows the CO J=3-2 distribution, especially in the region of 
the CO J=3-2 peak, e.g. in the overlap region and in the two nuclei. Conversion of the CO 
J=3-2 brightness to SCUBA equivalent flux and integration of the brightness contours in 
each map over their common area shows that line emission accounts for approximately 46% 
of the integrated 850 /xm flux in the region shown in Figure 9(a). This is somewhat greater 
than the contamination factor of 30% estimated by Haas et al. (2000) based on a best 
estimate in the absence of actual CO J=3-2 data. The most serious contamination occurs 
in the southern gas complex of the overlap region (Kl in the map of Haas et al. 2000), 
where about 55% of the flux is from CO J=3-2. This region is coincident with the high 
excitation region seen in Figure 4(b) indicated by the r 3 i ratio. 

Figure 9(b) shows the result after subtracting the contamination by CO J=3-2 
emission from each pixel of the original SCUBA image. This image shows the emission 
from the dust alone overlaid on the CO J=3-2 greyscale map. The two 850 /xm knots in 
the overlap region (identified by Haas et al. as Kl and K2) are more widely and clearly 
separated, with the northern knot K2 being slightly stronger than the southern knot Kl. 
Kl suffers more contamination from the strong CO J=3-2 emission since it appears at the 
same strength as K2 in the uncorrected map. The corrected 850 /xm flux ratio between Kl 
and K2 S(K1)/S(K2) is 0.78. This ratio is similar to that for the 20 cm radio continuum, 
but is different from the ratio seen in the MIR which is 2.6 (Mirabel et al. 1998). The 450 
/xm flux from Haas et al. (2000) is also much higher in Kl than in K2 (no detection in K2), 
suggesting that the dust is hotter in Kl. As mentioned above, the CO J=3-2/J=l-0 ratio 
is also higher in Kl than in K2. 

Another striking feature in the corrected 850 /xm map is the lack of dust emission 
concentrated in the two nuclei. The cold dust emission seems to be diffusely distributed in 
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the disks of NGC 4038. In NGC 4039, we see almost no cold dust at all at the NGC 4039 
nucleus. Instead, more dust is found along the spiral arm southwest of the nucleus. This 
dust complex is detected in both 450 /im and 850 /im, and thus is most likely to be real. 
However, it is only a 4<r detection, and more sensitive measurements are needed to confirm 
it. 

Color plates 1 and 2 show the overlay of the corrected SCUBA 850 /im contours 
on the HST full color and Ha images from Whitmore et al. (1999). A dust lane in the 
overlap region is clearly seen in the HST full color image (Color plate 1). The SCUBA 850 
/im emission coincides with the dust lane, and young clusters are mainly seen at the edge of 
the dust lane. Some clusters are heavily obscured by dust. These clusters appear to be faint 
and red in the color image, but are clearly seen in the Ha image in Color Plate 2. Since Ha 
emission requires the presence of O and B stars to ionize the gas and these stars last less 
than 10 Myr, the presence of Ha emission alone guarantees that the region contains clusters 
younger than 10 Myr. The size of the Ha bubble can provide further information on the 
age of the clusters. The Ha bubbles within the 850 /im peak Kl and K2 are very compact, 
which suggest that this region is relatively young, while the cluster complexes bordering Kl 
and K2 and those along the northeast side of the dust lane have bigger bubbles and are 
probably older. Apparently, there has not been enough time for the complexes to blow away 
all the dust yet. Whitmore et al. (1999) estimated that the age of these clusters is ~ 5 Myr. 

In the Ha image (Color plate 2), there appear to be more clusters near K2 than Kl, 
and yet Kl generates much more MIR and NIR emission than K2. Thus the clusters in K2 
may be so young that the star formation is still deeply cocooned and hence almost invisible 
at optical or even NIR and MIR wavelengths. The age of these clusters should be less than 
5 Myr. 

Figure 10 shows contours of the corrected SCUBA 850 /im (white) and 450 /im (black) 
emission overlaid on the Nobeyama 45m CO J=l-0 greyscale map. The differences 
between the dust emission and CO emission are remarkable. In the overlap region, the 
CO J=l peak is located between the two 850 /im peaks, which is similar to the CO 
J=3-2 distribution. The 850 /im emission in the two nuclei also seems to be offset from 
the CO J=l-0 peaks, especially for the dust concentration southwest of the NGC 4039 
nucleus where little CO J=l-0 emission is detected in the Nobeyama single dish or OVRO 
interferometer maps by Wilson et al. (2000). However, we do see weak extended features 
in the BIMA + NRAO 12m full synthesis CO J=l-0 map (Gruendl et al, in preparation) 
in the position coincident with the dust emission. The difference between the CO and dust 
850 /im continuum suggests that the molecular gas traced by CO is not closely associated 
with the cold dust emission. This is consistent with our conclusion from the excitation 
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analysis in which the 12 C0 J =1-0 emission is mainly from a warm and moderately opaque 
component. 

Table 5 shows the total dust mass in the system computed using the corrected 850 
/xm flux. The total flux is well fitted by a one-component dust model with an emissivity 
law A~ 2 (p = 2.0), dust temperature T=28 K and mass M dust = 1.7 x 1O 7 M . The new 
estimate of dust mass is only 20 percent of the value suggested by Haas et al. (2000) 
(~ 10 s M ), which is attributable to our larger estimate of the contamination by CO J=3-2 
emission. Combining this with the gas mass of 2 x 1O 9 M , we obtain a gas to dust mass 
ratio of about 100. 

To fit the SED for the overlap region, mainly the two sub-mm knots, we assume that 
the FIR fluxes at 60, 90, 120, 160 and 200 /im in this region are 35% of the total flux from 
Klass et al. (1997). This fraction is estimated by Haas et al. (2000) from the ISOPHOT 60 
fim and 100 fim maps. A one-component model can fit the SED reasonably well (using 3a 
for K2 as an upper limit), with dust temperature T=28 K and mass M dust = 6.0 x 10 6 M . 
This mass is a factor of 4.4 times less than that estimated by Haas et al. (2000). We have 
also tried to fit the SEDs of the individual knots with the assumption that FIR emission is 
distributed equally between Kl and K2. A dust temperature of 28 K and (3 =2.0 can fit K2 
well within the error range, but Kl is better fitted with a dust temperature of 35K, j3 = 1.0 
for A < 250yum and (3 = 2.0 for A > 250/im (Hildebrand 1983). The estimated dust mass is 
listed in Table 5, with K2 containing slightly more dust mass than Kl. 



9. Origin of the molecular gas in the overlap region 

First of all, referring to the whole region between the two nuclei as "overlap" may be 
misleading. In §5 we showed that the kinematic data at the CO J=l-0 and 3-2 transitions 
suggest that the actual overlap occurs near the OV-line DEC(J2000) = — 18°52 / 35". Only 
the region near this line could have gas from both disks or in a bridge between two disks. 
The majority of CO and dust emission appears to be from the off-nucleus gas concentration 
in the disk of NGC 4039. If the conventional X factor is used, this concentration contains 
~ 3.2 x 10 9 M Q of molecular gas, more than the total M H2 in the Milky Way. This region 
also contains more than 30% of the total MIR, FIR, sub-mm and radio continuum emission, 
which makes it outshine the two nuclei at all these wavelengths. This is the place that 
contains most of the effects related to galaxy collision and the induced star formation. We 
refer to it as the "southern overlap" region to indicate its relation to the actual overlap 
region. The origin of the molecular gas in this region is crucial for understanding the 
star formation activity and for predicting the further evolution of this on-going merger. 
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Our LVG analysis indicates that the molecular gas may be ~ 5 — 10 times smaller than 
that estimated using a conventional X factor. In the following we argue that the lower 
value provides a reasonable estimate of the molecular gas based on the probable dynamical 
history, the short lifetime of molecular cloud and the small age constrained by the young 
clusters found in this region. 

The dynamical models suggest that following the initial encounter, the two galaxies 
separated and then later re-engaged to form the current configuration {e.g. Toomre & 
Toomre 1972; Dubinski et al. 1996). The overlap region is not the gravitational center and 
the "overlap" is just a projection phenomenon. Within the dynamical time scale (10 7 - 
10 8 years), proper motion and rotation of the disks would cause a relative displacement 
of the overlapping gas from the two galaxies. Hence it is difficult to understand how this 
particular region outside the galactic nucleus could accumulate large amounts of molecular 
gas. If this is an effect of the tidal force, we should also see similar effects in the stellar 
or atomic gas components in the same region, which is not the case. On the other hand, 
even if the accumulation of gas were possible dynamically, such a large scale dynamical 
process in a disturbed interacting system would unavoidably trigger massive starbursts in 
the GMCs which would then quickly destroy the molecular clouds. Hence the accumulation 
process of the molecular gas would not last longer than the life time of molecular cloud, 
which is ~ 10 7 years. Indeed, a number of young clusters have been found in the overlap 
region and the age of the clusters was estimated to be ~ 5 Myr (Whitmore et al. 1999). 
This is consistent with the time scale we estimated if those clusters were formed during the 
process of molecular gas formation/accumulation. Taking the maximum sound speed of 
10-100 kms" 1 for the molecular clouds in the galactic disk, the gas can move only 0.1-1 
kpc in 10 7 years. Thus not much gas can be accumulated in one region during such a short 
period of time. 

Haas et al. (2000) explained the dust /gas in the overlap region with a "traffic jam" 
model, which suggested that initial inelastic cloud collisions would build a kernel into which 
the other clouds of the rotating disks are running, and making huge amounts of gas pile 
up together. This mechanism, however, can only build up to 10% of the gas mass, simply 
because the lifetime of the molecular clouds is less than 10 years which is only one-tenth 
of the typical rotational period for spirals. Hibbard et al. (2001) suggested that the overlap 
region represents material transferred from NGC 4038 to NGC 4039. Such a model can 
explain the blue-shifted velocity pattern in the northern overlap region. However, from 
our CO J =1-0 and 3-2 kinematic data, the average blueshifted velocity is 150 kms -1 and 
therefore the gas can move only about 1.5 kpc in 10 Myr. Therefore only a modest fraction 
of gas in the overlap region could have resulted from accretion from NGC 4038. 
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Consequently the SGMCs in the southern or northern overlap region and the two dust 
knots seen in SCUBA 850 /im most likely formed independently out of material originating 
nearby. If the gas distribution in the progenitor NGC 4039 is originally normal, the vicinity 
of the southern overlap should not have more than one-tenth of the total molecular gas of 
that galaxy, which is ~ 1O 9 M . The total HI gas in the Antennae galaxies is 6.3 x 1O 9 M , 
but 64% of this is associated with the tidal tails, and only one-third of the remaining HI 
flux is from NGC 4039 (Hibbard et al. 2001). Hence the local atomic gas originally in the 
south overlap region is also limited. If a large amount of H 2 is converted from HI, it is 
difficult to understand how to move a huge amount of atomic gas into the overlap region 
and convert it into H 2 in such an efficient way in a short time. 

In summary, the conventional X factor results in an extraordinary large amount of 
molecular gas outside the galactic nucleus which is difficult to explain with the existing 
theories of gas dynamics and molecular chemistry. On the other hand, if the X factor is 
actually 5-10 times smaller, the molecular gas mass would be ~ (3.2 — 16) x 1O 8 M which 
can be easily accounted for by the aforementioned mechanisms. Hence the strong CO 
emission in the southern overlap region is probably originated from the local molecular 
gas that is rendered more luminous by increased starburst activity. Our LVG analysis 
points to the possibility that large velocity dispersion and large filling factor of a diffuse gas 
component are associated with a low optical depth and high CO luminosity. The southern 
overlap is special only in that it is adjacent to the overlap region. Cloud collisions, shocks 
and induced star formation in this region can all contribute to the large velocity dispersion 
and create a diffuse component of molecular gas with excess CO J=l-0 emission. 

10. Summary 

We have obtained high quality, fully sampled, single dish maps of the Antennae 
galaxies NGC 4038/9 at the 12 CO J =1-0 and 3-2 transitions with an angular resolution of 
15" (1.5 kpc). This is so far the highest resolution single dish data at both the J=l-0 and 
3-2 transitions. 12 CO J=2-l data at the two nuclei as well as in part of the overlap region 
with 20" angular resolution have also been obtained. In addition, the 13 CO J=2-l and 
13 CO J=3-2 line emission are detected at selected points in the two nuclei and the overlap 
region. Our major findings are summarized below: 

I. The integrated 12 CO J=l-0 flux in the Nobeyama 45m single dish map is about 
twice that reported by Wilson et al. (2000) using the OVRO interferometer. The 12 CO 
J=l-0, 2-1, 3-2 emission all peak in an off-nucleus region adjacent to where the two disks 
overlap. This region accounts for approximately 42% of the total 12 CO J=l-0 flux of the 
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inner disk of the Antennae. 

A conventional X factor yields ~4x 10 9 M for the mass of molecular gas in the 
overlap region. Such large amounts of molecular gas are difficult to explain given the short 
lifetime of molecular clouds and transient overlap time scale. 

2. The ratio of the integrated intensities 12 CO J=3-2/J=l-0 (r 3 i) varies across the 
overlap region by a factor of 2. The two nuclei and southern overlap region have a higher 
r 31 ratio, which is coincident with the strong MIR and SCUBA 450 /im emission in this 
region. 

3. The ratio r 31 changes from less than 0.1 to as high as 0.9 over the velocity range 
from 1350 kms" 1 to 1620 kms" 1 near the OV-line DEC(J2000) = -18°52'35" which is the 
southern edge of the NGC 4038 disk overlapping with NGC 4039, indicating at least two 
velocity components with completely different excitation in this region. North of this line 
the excitation reverses between the high and low velocity components, i.e. the low velocity 
component becomes the high excitation one. The kinematics and excitation patterns 
suggest that the molecular gas north of this line is from NGC 4038 and south of this line is 
from NGC 4039. 

4. Both the R 2 i and i? 32 ratios are remarkably high in the overlap region. This is 
the first published case in which such high R 2 ± and _R 32 ratios are found outside a galactic 
nucleus. 

5. LVG models indicate that 12 CO and 13 CO emission comes from at least two different 
co-spatial components. The low density component has n(H 2 ) ~ 3 x 10 3 cm" 3 and is 
responsible for most of the 12 CO emission. This component must have a low to moderate 
optical depth (rio ^ 1), and/or a high [ 12 CO/ 13 CO] abundance ratio. The high density 
component has n(H2) ^ 1 x 10 4 cm" 3 and contributes more to the 13 CO emission. 

6. Assuming a Galactic [CO/H 2 ] abundance ratio of 10~ 4 , we find that the X factor for 
converting CO-to-H 2 is ~ 10 times smaller than the conventional value in the overlap and 
nuclear region of NGC 4038. Therefore, the molecular gas mass could be overestimated by 
a factor of up to 10 in the overlap region and in the two nuclear regions if the conventional 
X factor is used. However the mass of molecular gas derived by the LVG method used here 
is inversely proportional to the assumed fractional abundance of CO. 

7. The 12 CO J=3-2 line emission contributes significantly to the SCUBA 850 
fim continuum flux. For the Antennae galaxies, the ratio of the contribution from 12 CO 
J=3-2 emission to the SCUBA 850 fj,m flux is ~ 46% on average, and it varies by a factor 
of two across the system. After correcting for the 12 CO J=3-2 contamination, the dust 
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emission at 850 jum detected by SCUBA is consistent with the thermal emission from a 
single warm dust component with f3 — 2 and M^ust — 1-7 x 1O 7 M . This value is a factor of 
more than six times lower than that estimated by Haas et al. (2000) using a lower estimate 
for the contamination by 12 CO J=3-2 line emission. 

8. The corrected 850 /im continuum emission is not tightly correlated with the CO 
J =1-0 emission. This is perhaps related to the conclusion that the CO emission comes 
predominantly from a warm diffuse component which does not necessarily follow the 
distribution of the cold dust. The enhanced CO luminosity in the overlap region may be 
caused by this warm component heated by star formation. 
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Figure Captions 

Fig 1 The grid map of all the 12 C0 J=l-0 spectra. The temperature scale is 

T m b = —0.2. ..1.0 K, the velocity range is V LS r = 1240 — 1890 km s" 1 , and the resolution is 

AV chan = 10.5 km s" 1 . 

Fig 2 The grid map of all the 12 CO J=3-2 spectra. The temperature scale is 

T m b = —0.2. ..1.0 K, the velocity range is V LS r = 1161 — 1963 km s" 1 , and the resolution is 

AV chan = 10.8 km s" 1 . 

Fig 3 Comparison of the CO J=3-2 profiles (thick) with the J=l-0 profiles (thin). The 
former has been re-gridded and convolved to an angular resolution of 15" for comparison 
with the CO J=l-0 data. The temperature scale is T mb = —0.2. ..1.0 K, and the velocity 
range is V LS r = 1240 - 1890 km s _1 . 

Fig 4 (a) The 12 CO J=l-0 contour maps overlaid on the J=3-2 greyscale map (top). The 
contour levels are 16 K kms" 1 and increment by 16 K kms -1 . (b) The derived r 31 contours 
and greyscale map (bottom). The circle at the lower right corner indicates the resolution 
(15"). The contour levels are 0.10, 0.25, 0.4, 0.55, 0.70, 0.85. 

Fig 5 The 13 CO J=3-2 profiles overlaid on the re-observed 12 CO J=3-2 profiles at the 
center of the overlap region (the map center in Figures 1 and 2) and the NGC 4038 nucleus. 

Fig 6 The 13 CO J— 2-1 profiles overlaid on the 12 CO J— 2-1 profiles in the overlap region 
and the two nuclei. 

Fig 7 Channel maps of (a) 12 CO J=l-0 emission (top), and (b) 12 CO J=3-2 emission 
(bottom). For both maps, each velocity channel is 20 kms -1 . The offsets are relative to 
R.A.= 12 fe 01 m 54.8 s , DEC= -18°52'55" (J2000). The contour levels are 0.1, 0.2, 0.3. ..0.9 K 
in T m b units. A K-band NIR image is shown as the background of the last channel map. 
The star, triangle and cross symbols indicate respectively the CI, C2 and C3 components. 

Fig 8 The 12 CO J=3-2/l-0 contour maps overlaid on (a)the K-band NIR image (top) 
and (b) the ISO 15 /xm image (bottom). The circles at the lower right corner indicate the 
resolution which is 15". The contour levels are 0.10, 0.25, 0.4, 0.55, 0.70, 0.85. 

Fig 9 The SCUBA 850 /jm contours overlaid on the JCMT CO J=3-2 greyscale map (a) 
before (top) and (b) after (bottom) correcting for the CO J=3-2 contamination. The circle 
at the lower right corner indicates a 15" beamsize. The contour levels are 35, 58, 81, 104 
mJy/beam in (a) and 22, 35.5, 49, 62.5 mJy/beam in (b). 

Fig 10 The SCUBA 850 //m (white) and 450 fxm contours (black) overlaid on the Nobeyama 
CO J=l— map. The circles at the lower right corner indicate a 15" beamsize. The contour 
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levels are 22, 35.5, 49, 62.5 mJy/beam for the 850 fxm contours and 30, 50, 70, 90 percent 
of the peak for the 450 /jm contours. 

Fig 11 The corrected SCUBA 850 /im map (contours) overlaid on the HST full color image 
of the "Antennae" galaxies. The contour levels are 22, 35.5, 49, 62.5 mJy/beam. 

Fig 12 The corrected SCUBA 850 /j,m map (contours) overlaid on the HST H a image of 
the "Antennae" galaxies. The contour levels are 22, 35.5, 49, 62.5 mJy/beam. 
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Table 1: Integrated intensity ratios of 12 CO J=3-2/l-0 
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b The integrated intensity (T mb ) of CO J=l-0. 

c The integrated intensity (T m b) of CO J— 3— 2, which have been convolved from 14" to 15" and re-gridded to match the CO J — 1-0 data. 

The average integrated intensity ratios in a 15" beam. 
e The uncertainty in the r^i ratio. 
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Table 2: Observed line ratios a 



Position 


R21 b 


r 2 1 


R32 


»"31 


la r 32 


N4038 


13 ± 3 


1.01± 0.1 


19 ± 4 


0.91 ±0.1 


0.6± 0.2 


N4039 


16 ± 4 


1.19± 0.1 




0.73 ±0.08 




Overlap 












(0,0) 


26 ± 5 


1.04± 0.1 


25 ± 8 


0.78 ±0.09 


0.77±0.23 


(0,10) 


29 ± 5 


1.11± 0.1 




0.67 ±0.08 




(0,-10) 


18 ± 4 


1.06± 0.1 




0.87 ±0.1 





a The average line intensity ratios in a 20 beam. The intensities with a higher resolution have been convolved to 20 . 
b The quantities R21, r2i, -R32, T31 and 13 r32 arc defined in Section 6. 
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Table 3: One-component model fitting result 



Position 


n{H 2 ) 


T k 


X co /(dV/dR) 
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10 3 cm~ 3 
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43 
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0.018 


40 
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1- 10 
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2.2-3.1 


0.3-1.0 




40 


<2 


Overlap 


















(0,0) 


3.0 


36 


3.0 


2.5 


0.90 


0.051 


70 


2.2 




1- 6 
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1-10 


1.8-3.1 


0.3-1.0 




40-200 


<3 


(0,10) 


3.2 


33 


1.0 


0.98 


0.33 


0.061 


60 


1.2 




1- 5 


25-130 


1-10 


1.5-3.1 


0.2-1.0 




>40 


<2 


(0,-10) 


10 


25 


1 


3.1 


1.3 


0.030 


70 


0.9 




1- 20 


20-130 


1-10 


3.1-6.2 


1.0-1.5 




>40 


<2 



Note: for each position, the first row is the best fit parameters and the second row is the estimated 
range within the observed uncertainty. 
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Table 4: Two-component modeling for the overlap region 



\ rz COi 


LVG parameters a 

-■-kin) 11 \ ±± 2Ji lv i J 


Line ratios b 

7 o 1 l o 1 Pi o 1 Pi 

A 21-) i oli iL oz 


N CO /AV c 
ctd,~^ 1 km s ^ 


v • d 

A mm 


70 


(A) 36, 3.2xl0 3 , 3xl0" 6 , 0.051 
(B) 90, l.OxlO 5 , 3xl0~ 6 , 0.00033 


1.06, 0.80, 23, 28 
(1.04, 0.78, 26, 25) 


1.9 xlO 15 


0.9 


70 


(A) 120, l.OxlO 3 , lxlO" 5 , 0.044 

(B) 90, l.OxlO 5 , 3xl0" 6 , 0.0005 


1.07, 0.78, 24, 25 
(1.04, 0.78, 26, 25) 


1.8 xlO 15 


0.6 


70 


(A) 120, l.OxlO 3 , lxlO" 5 , 0.044 

(B) 42, 3.2xl0 4 , lxlO" 6 , 0.0022 


1.07, 0.79, 23, 29 
(1.04, 0.78, 26, 25) 


1.6 xlO 15 


1.1 


40 


(A) 120, l.OxlO 3 , 3X10" 6 , 0.10 
(B) 220, l.OxlO 5 , 3xl0" 6 , 0.00060 


1.13, 0.70, 23, 27 
(1.04, 0.78, 26, 25) 


1.5 xlO 15 


1.6 



a LVG parameters for component A (first row) and B (second row). / is the filling factor of each component. The units are the same ones used in 
Table 3. 

b The line ratios estimated for the combination of two components. The ones in the parentheses are the observed values. 

c Nco is the total column density combining the components using the formula: 

Nco/AV - [n{H 2 )X C o/(dV/dR) f] A + [n(H 2 ) Xco / (dV/dR) f] B 

d The Xmin parameter characterizes the goodness of the fit (see Table 3) 
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Table 5: Comparison of SCUBA 850 fim flux and derived dust mass before and after correcting the CO J=3— 2 contamination 



Position 


S'450 a 

Jy 


Sgso (Haas) 6 
mjy 


Ssso (revised) c 
m.Ty 


M dust (Haas) " 
10 7 M o 


Mdust (revised) ° 
10 7 M 


NGC4038 


1.69 


250 


166.8 


2.8 


0.5 


NGC4039 


1.06 


150 


< 78 


1.7 


<0.3 


Kl 


1.33 


190 


83.8 


>2.2 


0.25 


K2 


< 0.90 


170 


106 


>3.8 


0.30 



a The SCUBA 450 flux from Tabic 1 of Haas et al. (2000). 
b The SCUBA 850 flux from Tabic 1 of Haas et al. (2000). 

c The revised 850 /im continuum flux after correcting for CO J— 3-2 contamination. See text for details. 

The dust mass estimated by Haas et al. (2000). 
c The dust mass estimated from the real 850 /xm continuum flux. The dust model was the same as that of Haas et al. (2000). See text for details. 
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